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1. Introduction 

It is known that extremal configurations in many important problems in classical 
complex analysis exhibit complicated fractal structure. This makes such problems 
extremely difficult. The classical example is the coefficient problem for the class 
of bounded univalent functions. Let = z + a2z'^ + 03^'^ + . . . be a bounded 
univalent map in the unit disc. One can ask what are the maximal possible values of 
coefficients a„, especially when n tends to infinity. We define 7 as the best possible 
constant such that |a„| decays as n''~^. In Carleson and Jones showed that this 
problem is related to another classical problem about the growth rate of the length 
of Greens' lines. In particular, they showed that the extremal configurations for 
both of this problems should be of a fractal nature. 

During the last decade it became clear that the right language for these problems, 
as well as many other classical problems, is the maltifractal analysis. It turned out 
that all these problems could be reduced to the problem of finding the maximal 
value of the integral means spectrum. 

In the recent paper [2] S. Smirnov and the author introduced and studied a new 
class of random fractals, the so-called random conformal snowflakes. In particular, 
they proved the fractal approximation for this class, which means that one can find 
conformal snowflakes with spectra arbitrary close to the maximal possible spectrum. 
In this paper we report on our search for snowffakes with large spectrum. 

The paper organized as follows: in the introduction we give some basic informa- 
tion about integral means spectrum, define random conformal snowflakes, and state 
the main facts about spectrum of the snowflakes. In the Section[2]we give numerical 
estimates of the spectra of several snowflakes for different values of parameters. In 
the last Section we give rigorous lower bound for the spectrum at t = 1. 

1.1. Integral means spectrum. Here we brieffy sketch the necessary definitions 
and the background. For the detailed discussion of the subject we recommend 
surveys [3 [HI [1] and books p",?]. 

Let ri G C be a simply connected domain in the complex sphere which contains 
infinity. By the Riemann uniformization theorem, there is a conformal map from 
the complement of the unit disc onto f2 such that 0(cxd) = 00. The integral means 
spectrum is defined as 

log fn^'' \(b'(re">)\'d9 
r^i+ |log(r-l)| 
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The universal integral means spectrum is 

B{t) = sup/3(t), 

where supremum is over all simply connected domains fl. 

The universal spectrum plays the central role in the Geometric Function Theory. 
A lot of major problems in the field can be stated in terms of B{t). In particular, 
Brennan's conjecture [3] asserting that IV''!^"'^ is integrable for any e > and any 
conformal t/; : ^ D is equivalent to B{—2) = 1. The Carleson- Jones conjecture, 
7 = 1/4, is equivalent to B{1) = 1/4. In 1996 Kraetzer formulated the ultimate 
conjecture on integral means: 

B{t) = t^A, \t\<2, 
B{t) = |t|-l, |i|>2. 

This conjecture is based on the above mentioned conjectures, computer experiments 
and the believe that the spectrum should be relatively simple function. Kraetzer 
performed computer experiments with quadratic Julia sets. The results were within 
5 percents from i^/4. We should point out that the experiments were completely 
non-rigorous, and, as far as we know, the error can not be estimated using available 
techniques. 

In the present paper we give a computer assisted lower bound on the universal 
spectrum which are within 4— 7 percents from t^/4 but they are "semi-rigorous" . By 
this we mean that we use computer only for numerical integration. For numerical 
integration of an explicit function one can write estimates of the error term and 
this will give a rigorous estimate. We give an estimate of the error term only for 
t = I (this case is of special interest and estimates are a bit simpler). We prove that 
B(l) > 0.23 which is a significant improvement over previously known B[\) > 0.17. 

Essentially the same argument as for the case t = 1 could be used to estimate 
the error for other values of t, although we are not doing it here. 

1.2. Conformal snowflake. We denote by S' the class of univalent maps cj) from 
the complement of the unit disc D_ into itself such that 0(oo) = oo and (p'ioo) e M. 
Let Kn<j) be the Koebe transform of S S': 

(/^„</))(z) = (i^„0)(z) = 'VH^. 

We denote the conjugation by rotation by (peiz) — e*^(/)(ze^*^). 

To construct a random snowflake we need two components: a building block 
e E' and an integer k > 2. The n-th approximation to the snowflake is 

fn{z) = /„-i(iffc.0e„(^)) = <l>e„{4i\- ■ ■ '^ef (^'") • ■ • ), 
where 9j are independent variables uniformly distributed between and 27r. The 
conformal snowflake / is the limit of /„. For the proof that conformal snowflakes are 
well defined see [2]. 5 = C \ /(©_) and g = f^^ are also referred to as snowflakes. 

It turns out that the spectrum of a snowflake for a particular value of t is closely 
related to the spectrum of an integral operator P whose kernel is defined in terms 
of fc, t, and (j). This operator is defined as 

Let R be any radius such that Dji C tp'^{Dji) where Dji = {z : 1 < |z| < R} and 
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Figure 1 . The image of a small boundary arc under /3 with three 
Green's lines. 



One can show (see [2] for the proof) that if A is the maximal eigenvalue of P 
then 



(2) 

in particular 
(3) 



>logA/logfc, 



/3W> mm^log(^)/logfc 



for any positive test function i>. 

1.3. Random fractals. One of the main problems in the computation of the inte- 
gral means spectrum (or other niultifractal spectra) is the fact that the derivative 
of a Riemann map for a fractal domain depends on the argument in a very non reg- 
ular way: (p' is a "fractal" object in itself. The solution to this problem is to study 
random fractals for which the boundary behavior of is statistically independent 
of the argument. In this case it is natural to consider the average integral means 
spectrum: 



m = sup /(r-l) 



N/3-1 



= inf |/3 



(r-1) 



/3-1 





27r 



^ [\fire'^)\*] dOdr = oo 
IE [|/'(re'^)r] dedr < oo 



The average spectrum does not have to be related to the spectra of a particular 
realization. We want to point out that even if (jj has the same spectrum a.s. it does 
not guarantee that f}{t) is equal to the a.s. value of f3{t). Moreover, it can happen 
that (3 is not a spectrum of any particular domain. On the other hand, Makarov's 
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fractal approximation ([5]) implies that P{t) < B{t). Thus lower bomid for any (3 is 
a lower bound for B. In the remainder of the paper we will work only with average 
integral means spectrum. We refer to it as "the spectrum" and denote it by [3. 

2. Numerical estimates of the spectrum 

In this section we give a numerical support of Kraetzer's conjecture. To achieve 
this we estimate the spectrum of snowflakes with a straight slit building blocks 
(see Figure 1). We use the following scheme: first we compute the disctretized 
operator P and find its main eigenvalue and eigenvector. We approximate the 
main eigenvector by a nice function v and use it as a test function in 

First we define the basic slit function 



(4) (l){z,l)^(l)i{z)^ fl2 



V/if(z.s) + /V(4fc + 4) 



where s is a constant close to 1, fj,i and fi2 are the Mobius transformation that 
maps D_ onto the right half plane and its inverse: 

z-1 



z + V 
z + l 



z- 1 

We also need the inverse function 

(5) v(^,0 = V'iW = 0(^,0"'- 

The image (j)i{IS>-) is D_ with a horizontal slit originating at 1. The length of the 
slit is I. The derivative of a slit map has singularities at points mapped to 1. But 
for s > 1 these singularities are not in D_. 

Given k, I (and s) we can numerically find the critical radius R and define 
the operator P. For the first step, we use the following discrete operator which 
approximates P. Choose sufficiently large N and M. Let r„ — 1 + (i? — l)n/N and 
9m = 27rm/Af . P is approximated by an A'^ x matrix with elements 

|0'(ry^'«-)|* 



Pn,n' — ^ ^ ^ 



-t{k-l)/k_ 



where summation is over all indexes m such that r„' is the closest point to |(/)(ry^e'^™)|. 
This defines the discretized operator P^. Let Xn and Vn be the main eigenvalue 
and the corresponding eigenvector. 

It is highly plausible that the operators Pn converge to P and the A n converge 
to the main eigenvalue of P, but it seems to be hard to prove. Moreover, the rate 
of convergence is difficult to estimate. But this crude approximation gives us the 
fast test on whether the pair (jj and k defines a snowflake with large spectrum or 
not. 

For each value of t we use this method to search for the optimal values of k and 
I. This gives us snowflakes that a supposed to have large spectrum. To estimate 
the spectrum we use the following trick. We approximate the eigenvector Vn by a 
simple function (or just linearly interpolate it). The result is used as a test function 
v. Than we use numerical integration to calculate Piy. By ^ it gives us the lower 
bound of p. For general values of t we do it without estimates of the error term, we 
just use Euler quadrature formula doubling number of nodes until the difference is 
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t 


k 


I 


log I, A 


0it) 


Kraetzer 


tV4 


-2.0 


34 


1 


1.262 






1.0 


-1.8 


34 


1 


1.068 






0.81 


-1.6 


34 


1 


0.8761 






0.64 


-1.4 


34 


1 


0.6879 




0.476 


0.49 


-1.2 


34 


1 


0.5059 




0.340 


0.36 


-1.0 


34 


1 


0.3354 




0.231 


0.25 


-0.8 


34 


1 


0.1865 




0.149 


0.16 


-0.6 


24 


21 


0.0848 


0.0710 


0.085 


0.09 


-0.4 


20 


25 


0.0377 


0.0352 


0.037 


0.04 


-0.2 


31 


44 


0.0093 


0.0083 


0.0095 


0.01 


0.2 


5 


7 


0.0091 


0.00897 


0.0094 


0.01 


0.4 


11 


30 


0.0376 


0.03767 


0.037 


0.04 


0.6 


14 


68 


0.0851 


0.08442 


0.086 


0.09 


0.8 


12 


67 


0.1514 


0.1511 


0.154 


0.16 


1.0 


13 


73 


0.2362 


0.2340 


0.242 


0.25 


1.2 


10 


67 


0.3425 


0.3350 


0.346 


0.36 


1.4 


8 


55 


0.4680 


0.4586 


0.476 


0.49 


1.6 


6 


39 


0.6137 


0.6091 




0.64 


1.8 


6 


39 


0.7790 


0.7713 




0.81 


2.0 


4 


21 


0.9548 


0.9296 




1.0 



Table 1 . Spectra of nearly optimal snowflakes for different values 
of t. Parameters k and I define the snowflake (s=l), A is the main 
eigenvalue of Pn (computed with N 2000 and M = 1000), (3(t) 
is the lower bound given by ([3]). In the last two columns we give 
lower bounds obtained by Kraetzer and values of the conjectured 
universal spectrum. 



small. We also can argue that since functions seems to be quite smooth (there is 
no oscillation), the precision of numerical integration should be much better than 
what follows from the standard estimates of the error term. The results of our 
computations are given in the Table 1. 

We want to make a few comments about negative values of t. The derivative of 
(j) has a zero of the first order at z = 1, this means that for p = 1 the integral in the 
definition of P diverges at r = 1. This is the reason why we can not find a good 
test function and the estimated A is probably far from the correct one for t < —0.8. 
The problem with A could be solved by increasing N and M . For example if we 
set N = 3000 and M = 2000 then for t = -0.6, . . . , -1.2 the logarithm of A will be 
equal to 0.0847, 0.1579, 0.2926, and 0.4509 correspondingly. This means that the 
value for t — —0.6 is probably correct even for smaller values of N and M and the 
value for t = —1.2 is still far from its true value. 

3. Spectrum at t = 1 

In this section we estimate the error term in the numerical integration of Pv. 
This gives a rigorous proof that B{1) > 0.23 which is a significant improvement 
over previously known B{1) > 0.17. We also would like to point out that this 
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computation is not specific for t = 1 (though it is a bit easier in this case) and can 
be done for other values as weU. 

For t = 1 we study the snowflake with / = 73 and fc = 13. The Figure [T] shows 
the image of a small arc under /a and three Green's lines. This snowflake has a 
spectrum close to 0.234. Unfortunately there is a singularity in the kernel of P 
(it is integrable, so one can still estimate the error term). This singularity is the 
reason why we introduced the parameter s. If we choose s close to 1 then snowflake 
is close to the snowflake with s = 1, but derivative of is continuous up to the 
boundary. In our case we use s — 1.002. 

First we have to find the critical radius R such that Dj^ C ijj^{Di^). By symmetry 
of (j), the critical radius is the only positive solution of 

■0''(x) = X. 

This equation can not be solved explicitly, but we can solve it numerically (we don't 
care about error term since we can use any upper bound for R). The approximate 
value of R is 76.1568. To be on the safe side we fix i? = 76.2. 

To find the test function we approximate P by Pjv where N = 1000 and M — 500. 
The logarithm of the first eigenvalue is 0.2321 (it is 0.23492 if we take s = 1). We 
scale the coordinates of the main eigenvector of Pn from [1, 1000] to the interval 
[1, R] and approximate by a rational function v which we will use as a test function: 

uix) = (7.1479 + 8.9280X - 0.07765^^ + 1.733 x lO^^x^- 
2.0598 X IQ-^a;'' + 9.5353 x 10~^a;''^)/(2.7154 + 13.2845a;). 

To estimate (3{1) we have to integrate :^(|0|)|(/)'|/|(/)|. It is easy to see that the 
main contribution to the derivative is given by a factor |0'|/|0|. Assume for a while 
that s = 1. The fraction |0'(z)/0(z)| can be written as 

(6) l^-'l 



l^l\/|(2-2l)(2-2:2)' 

where 



. Z^^^^H^^ . -0.894756 - 0.446556. 
5625 



-5033 + 292iV74 

Z2 = — w -0.894756 + 0.446556?. 

5625 

Singular points zi and Z2 are mapped to 1 and has a square root type singu- 
larity at these points. They will play essential role in all further calculations. We 
introduce notation zi — x + iy and Z2 = x — iy, for z we will use polar coordinates 

iO 

z — re . 

We compute the integral oi f — |i^((/))0'/0| using the Euler quadrature formula 
based on the trapezoid quadrature formula 

/■27r 

/ f{x)dx « S:{f) = S,{f) - 5]72fee'' (/(^'=-i)(2^) - f^^'-'Ho)) , 

k=i 

where S({f) is a trapezoid quadrature formula with step e and jk = Bk/k\ were Bk 
is the Bernoulli number. The error term in the Euler formula is 

(7) -72„max/(2")e2n27r. 
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In our case function / is periodic and terms with higher derivatives vanish. This 
means that we can use ([7]) for any n as an estimate of the error in the trapezoid 
quadrature formula. 

Function has two singular points: zi and Z2. Derivative of <j) blows up near 
these points. This is why we introduce scaling factor s. We can write a power series 
of <f) near zi (near Z2 situation is the same by the symmetry) 

This means that for s > 1 derivative can be estimated by 



Co 



1) 



-k+l/2 



+ \ck+i\{s-l) 



-fc+3/2 



The series converges in a disc of a fixed radius (radius is |zi + 1|), so its tail can be 
estimated by a sum of a geometric progression. Writing these power series explicitly 
we find (for s = 1.002) 

|0'| < 55, |(/)"| < 11800, 10(3^1 < 8.69 X 10^ 

IcZ-^"*)] < 1.08 X 10^", < 1.90 X 10", l./)^^)! < 4.25 X 10l^ 

|(/.(^)| < 1.17 X 10^0. 

The maximal values for first six derivatives of v are 

|i^'|<0.28, \v"\<QAb, < 1.12, 

< 3.69, < 15.3, < 76.2. 

The derivative de\(f)\ can be estimated by r\(t)'\. We can write sixth derivative of 
j^(|(/)|)|0'|/|0| as a rational function of partial derivatives of |(/)|, |(/)'|, and v. Than 
we apply triangle inequality and plug in the above estimates. Finally we have 



d 



'^(101)10' 



< 1.65 X 10^1. 



Plugging in the value e — 7r/5000 and the estimate on the sixth derivative into ([7]) 
we find that error term in this case is less than 0.0034. 

Next we have to estimate modulus of continuity with respect to r. First we 
calculate 

dr\z - (a + 6i)P = 2r- 2(acos6' + 6sin6'). 
Applying this formula several times we find 



k-ll 



101 J \r^z - zi\\z - Z2\ 

r — cos 9 r ~ X cos 9 ~ y sin C 



< dr 



\z-l\ 



MS 



2\z~zi\^S 



z-1 



\\/\z - Zl\\z - Z2 

r — X cos 9 + y sin " 



2\z^Z2\^S 



\z-l\ 



where S — ^\z — zi\\z — Z2\- Factoring out 

1 



2|z-l|-|z-Zi|5/2|2-Z2p/2 



we get 
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2(r - cos6')|z - zip|z - ~ {r ~ xcosO ~ ysm9)\z - - 

— (r — xcosd + y sin0)|z — — Zip 
= -2(r2 - l)(2cos2 6'r(a; - 1) + cos9{r'^ + l)(x - 1) + 2ry2). 
This is a quadratic function with respect to cos0. Taking values of x and y into 



account we can write it as 



cos 6 + cos ( 



1\ 1 592 



2 5625 



This quadratic function has two real roots. Their average is — (r + l/r)/2 < — 1, 
hence one root is definitely less than —1. The product of roots is a small negative 
number, which means that the second root is positive and less than 1. Simple 
calculation shows that this root decreases as r grows. This means that the corre- 
sponding value of 6 increases. Hence it attains its maximal value at r = 1.4 and the 
maximal value is at most 1.48. This gives us that the radial derivative of |^'|/|0| 
can be positive only on the arc 9 € [—1.48, 1.48]. By subharmonicity it attains the 
maximal on the boundary of {z | 1 < r < 1.4, —1.48 < 9 < 1.48}. It is not very 
difficult to check that maximum is at z = 1.4 and it is equal to 0.36. 
Let 



0'(r 



l/k i6 



d9 
2^' 



J —TT 

The derivative is 

By the symmetry the first integral is zero. In the second integral 

can be positive only when 9 e [—1.48, 1.48] and even in this case it is bounded by 
0.36/(ri-i/'=/c27r). Hence 

I'{r) < 2 • 1.48 • 0.36/(r^"^/'''fc27r) < 0.0131r^/'=~^ 

If we compute values I{ri) and I(r2) (with precision 0.0034) then the mininmm of 
P{v)/v on [ri,r2] is at least 

(8) rj/''(min{/(ri),/(r2)} - 0.0034 - 0.0131(r2 - riy/^-^)/v{ri). 

We take 3000 equidistributed points on [1, R] and compute I{r) at these points. 
Applying the error estimate ([8|) we find a rigorous estimate from below of Pv jv. 
The minimum of Pvjv is at least 1.8079 which means that 

/3(1) > 0.2308. 
References 

[1] D. Beliaev and S. Smirnov. Harmonic measure on fractal sets. Proceedings of the 4th European 

congress of mathematics, 2005. 
[2] D. Beliaev and S. Smirnov. Random conformal snowflakes. Preprint, 2006. 
[3] J. E. Brennan. The integrability of the derivative in conformal mapping. J. London Math. Soc. 

(2), 18(2):261-272, 1978. 
[4] L. Carleson and P. W. Jones. On coefficient problems for univalent functions and conformal 

dimension. Duke Math. J., 66(2):169-206, 1992. 



INTEGRAL MEANS SPECTRUM OF RANDOM CONFORMAL SNOWFLAKES 



9 



[5] N. G. Makarov. Fine structure of harmonic measure. St. Petersburg Math. -J., 10(2);217— 268, 
1999. 

[6] C. Pommerenke. Univalent functions. Vandenhoeck & Ruprecht, Gottingen, 1975. 

[7] C. Pommerenke. Boundary behaviour of conformal maps, volume 299 of Grundlehren 

der Mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences]. 

Springer- Verlag, Berlin, 1992. 
[8] 0. Pommerenke. The integral means spectrum of univalent functions. Zap. Nauchn. Sem. S.- 

Peterburg. Otdel. Mat. Inst. Steklov. (POMI), 237(Anal. Teor. Chisel i Teor. Punkts. 14) : 119- 

128, 229, 1997. 



